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A negatively charged nitrogen vacancy (NV) center in diamond has been recognized as a good 
solid-state qubit. A system consisting of the electronic spin of the NV center and hyperfine-coupled 
nitrogen and additionally nearby carbon nuclear spins can form a quantum register of several qubits 
for quantum information processing or as a node in a quantum repeater. Several impressive exper¬ 
iments on the hybrid electron and nuclear spin register have been reported, but fidelities achieved 
so far are not yet at oor below the thresholds required for fault-tolerant quantum computation 
(FTQC). Using quantum optimal control theory based on the Krotov method, we show here that 
fast and high-fidelity single-qubit and two-qubit gates in the universal quantum gate set for FTQC, 
taking into account the effects of the leakage state, nearby noise qubits and distant bath spins, can 
be achieved with errors less than those required by the threshold theorem of FTQC. 
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I. INTRODUCTION 

Nitrogen vacancy (NV) centers in diamond have many 
remarkable properties. For example, the spins of NV cen¬ 
ters have relatively long relaxation and coherence time 
(even at room temperature) M, and the electron spin 
triplet ground state can be initialized, manipulated and 
read out with microwaves and lasers. These ex¬ 

ceptional properties make the NV center(s) a promising 
system for sensitive magnetic field and weak signal sens¬ 
ing [lol - l^ . bio-marking tracking (2ll - l^ and quantum 
information processing [3, [25l - l50j| . 

Quantum gate operations in a quantum register of in¬ 
dividual electron spin or/and nearby individual nuclear 
spins associated with a NV center in diamond have been 
demonstrated experimentally li^. However, 

the gate fidelities in theses experiments or studies are 
limited to certain values because the pulses sequences to 
perform the gates even in the ideal unitary case are not 
optimally designed and come with some sorts of approx¬ 
imation. 

There have been schemes proposed for protecting 
quantum gates of NV center spins from decoherence, 
based on dynamical decoupling (DD) protocols or/and 
dynamically corrected gate (DCG) [lo, SI, Slj Isil - ls^ . 
Experimental realizations of noise-resilient or decoher¬ 
ence protected quantum gates on NV centers have been 
reported [dlj. So far, only the single-qubit gates are 
shown to be of high-fidelity. For example, the fidelity 
of a dynamical-decoupling-protected X gate is shown to 
be about 0.985 for a gate duration of 35.5 fis [ij and the 
fidelity of a SUPCODE 7r/2 gate is about 0.9961 with a 
gate time of 5.063 fis Sl|- However, these gate times are 
much longer than those of unprotected gates, and how 
to practically implement the different protected gates for 
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different qubits in parallel in a many-qubit register is not 
clear. 

An alternative approach to realize high-fidelity quan¬ 
tum gate is through the quantum optimal control (QOG) 
|54 - l71| . A recent study [5^ investigated the theoretically 
achievable fidelities when coherently controlling an effec¬ 
tive three-qubit system consisting of a negatively charged 
(i^NV) center in diamond with an additional nearby 
carbon nuclear spin. The results in this study in¬ 
dicates that by using square and two frequency com¬ 
ponent radio and microwave frequency pulses, the best 
single-qubit gate fidelity is less than 98% and the multi¬ 
qubit gate fidelities are somewhat lower than that. It 
was thus suggested that to reach the fidelity threshold(s) 
predicted by current models of fault-tolerant quantum 
computation (FTQG) [73 - l75l |. going beyond the square- 
pulse paradigm and using pulse-shaping techniques like 
optimal control is required [H^. The QOG theory has 
been applied to NV center based quantum information 
processing (dTl - li^ . Reference!!^ considered only the elec¬ 
tron spin and designed fast single-qubit gates using the 
chopped random basis quantum optimization algorithm 
without resorting to the standard rotating-wave approx¬ 
imation condition. References 113 considered a system 
of a NV center’s electron spin and nitrogen ^“^N nuclear 
spins as well as coupled carbon nuclear spins, form¬ 
ing a small quantum register, and used the gradient as¬ 
cent pulse engineering (GRAPE) optimization algorithm 
[591461] to perform phase-flip quantum error correction 
on three qubits of one ^"^N nuclear spin and two nu¬ 
clear spins. Reference performed quantum gates and 
generated entangled states for two proximal NV centers 
in diamond using the GRAPE optimization algorithm. 
With the help of QOG, unwanted off-resonance transi¬ 
tions or crosstalk, and unwanted dipolar couplings be¬ 
tween the spins of the two proximal NV centers were 
significantly suppressed. However, in these QOG stud¬ 
ies and experiments [43 - 14^ . the maximum hyperfine in- 
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teraction strength between the NV electron spin and ei¬ 
ther the nitrogen nuclear spin or carbon nuclear spin is 
only about a few MHz. This is different from the case 
studied in Ref. where the hyperfine interaction of the 
15NV- is about 3MHz, while the hyperfine interaction 
with the nearest neighbor carbon can be larger than 
lOOMHz. A large hyperfine interaction potentially leads 
to fast quantum gate operations in the hybrid spin reg¬ 
ister. It is, however, this large hyperfine interaction that 
limits the maximum fidelity that can be achieved in the 
square-pulse paradigm [s^. Furthermore, these studies 
[471 - 14^ do not take the decoherence effect from the sur¬ 
rounding distant bath spins into the optimization consid¬ 
eration when constructing their QOC gates. 

In this paper, we present a detailed QOC study based 
on the Krotov optimization method I62 h^ m- 

ItH . [ 7 ^ for single-qubit and two-qubit gates of a hybrid 
electron and nuclear spin register of a NV center in di¬ 
amond taking into account the effects of leakage state, 
nearby noise spins and a bath of distant nitrogen spins 
or/and nuclear spins randomly distributed in the 
diamond lattice. In our model, a nuclear spin of 
which is in the first coordination shell (nearest neighbor) 
around the center and has strong hyperfine inter¬ 

action with the NV center electron spin is considered (see 
Fig.mfor a schematic illustration). The Krotov optimiza¬ 
tion method we employ has several ap pea ling advantages 
over the gradient methods [ssl - fs^ lo 1^ 1^: (a) mono¬ 
tonic increase of the objective with iteration number, (b) 
no requirement for a line search, and (c) macrosteps at 
each iteration. Quantum gates constructed via our QOC 
scheme with experimentally available or realistic param¬ 
eters are all with very fast speed and very high fidelity. 
Setting the gate operation times for our single-qubit X 
gate and Z gate performed on the electron spin to be 
10 ns, we obtain corresponding gate infidelities or errors 
to be 3.9 X 10“^ and 6.0 x 10“^, respectively. The two- 
qubit controlled-Not (CNOT) gate performed on the NV 
center electron spin and a proximal nuclear spin can 
be operated within 50 ns with an infidelity or error of 
about 4.3 X 10“"^ even in the presence of a host ^^N noise 
nuclear spin and an additional spin bath (environment) 
with a wide range of decoherence parameters. 

This paper is organized as follows. We briefly describe 
the model Hamiltonian of the NV-center-based hybrid 
spin register we consider in Sec. IH] To incorporate the 
effect of distant nuclear spins which form a spin bath or 
environment on the dynamics of the system qubits, we 
use the open-system master-equation approach and the 
description of this approach is presented in Sec. nn In 
Sec.IVl we define the infidelity or error function to mea¬ 
sure how well the gate operations of our system qubits 
deviate the ideal target gates in the presence of nearby 
noise qubits (spins) and a spin bath. The QQC algorithm 
based on the Krotov method is also briefly described here. 
In Sec. El we explore the application of the QQC for the 
implementations of Z-gate, A-gate and also the gates in 
the universal discrete quantum gate set for FTQC for the 
hybrid spin register. Comparison to the traditional ap¬ 




FIG. 1. (Color online) Schematic illustration of a negatively 
charged ^®NV“ center (V denoting a vacancy site) with one 
nearest neighbor atom, (a) For single-qubit operation, the 
electron spin of the NV center (in blue) is the system qubit 
while the ^®N and nuclear spins are regarded as noise 
qubits (in red), (b) For two qubit operations, the electron spin 
and the nuclear spin are regarded as the system qubits (in 
blue), while the nuclear spin is treated as a noise qubit (in 
red). The shaded area denotes additional distant bath spins. 


proach of implementing quantum gates and the effect of 
spin bath on the QOC gate operations are also presented. 
Finally, a conclusion is given in Sec. I VII 


II. MODEL HAMILTONIAN 

We consider a negatively charged NV center associ¬ 
ated with a ^®N nucleus (i.e., ^®NV“) in diamond . 
The electronic structure of the NV center has a spin- 
triplet ground state S' = 1 with a zero-held splitting 
A = 2.87 • 27rGHz between the mg = 0 and nis = ±1 
levels. Note that the quantization axis of this splitting 
is along the symmetry axis of the NV center, which we 
take as the z axis. The ^®N carries a nuclear spin I = 
Also, we consider a ^^C atom occupies one of the nearest 
position around the NV center, and other nuclear spins 
further away are regarded as a spin bath. A schematic 
structure of our system is shown in FigHJ Applying a 
static magnetic held B along the z axis splits the levels 
nis = — 1 and iris = 4-1. In our study, we consider all the 
three electron spin levels, i.e., the iris = 0 and nis = ±I 
levels, choose the spin levels nig = 0 and nig = —1 to be 
the two computational states of our electron spin qubit 
and treat the nig =4-1 state as an ancilla or a leak¬ 
age state. The NV center electron spin is coupled to the 
proximal ^®N and ^^C nuclear spins and an additional 
spin bath. 

It has been shown that the coupling of a spin bath of 
distant nuclear spins to a NV center electron spin can 
be modeled through classical magnetic held noise that 
causes decoherence by imprinting a random phase on the 
NV center electron spin [^. This semiclassical noise 
model emerges as the weak-coupling limit of quantum- 
mechanical entanglement-induced decoherence. The con¬ 
dition for this to be valid is when the electron-nuclear 
spin couplings quantihed by the magnetic helds Hnv.k 
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of the NV center electron spin at the sites of the nu¬ 
clear spins are much smaller than the externally applied 
magnetic field (i.e.,i?NV.n B^) [t^- In our study, 
we apply a relatively strong background magnetic field 
Bz to the spin register. We thus treat the coupling of 
the spin bath to the NV center electron spin to be ap¬ 
proximated as a classical random field B{t) which may 
depend on time acting on the z-component of the NV 
center electron spin, a pure dephasing model [zi- 
113. The couplings of the nearby or neighboring nuclear 
spins to the NV center electron spin are, however, treated 
quantum-mechanically. 

The total Hamiltonian of the system we consider can 
be written as follows. 

H = Ho + H,^ + H,y + H,b ( 1 ) 

Hq = AS'^ — 'yeBzSz — 'ycBzIcz — InBzInz 
+ AfSzIcz + Af + Sylcy) 

+ AfSziNz + Af {SJn. + SylNy) , (2) 

Bcxif) — Bx(t) ( 'JeSx 'yC^Cx 5 (3) 

Hcyit) = By{t) (—7e<S'y — y/clcy ~ y/NlNy) t (4) 
HeB = SzBf, (5) 

where Si is the spin-1 operators for the NV center elec¬ 
tron spin and Ici and Ini are the spin-1/2 operators 
for the and ^®N nuclear spins, respectively. 7 e = 
-2.8-27rMHz-G-i, 7 c = 0.00107 • 27r MHz • G"! and 
7 jv = — 0.43 • 27rkHz • G“^ are the gyromagnetic ratios 
for the electron, ^^G and ^^N spins, respectively [b^ . 

= Af = 127 ■ 27rMHz is the hyperfine coupling be¬ 
tween the electron and the nuclear spin [i^l, and 
Af = Af = 3.03 • 27rMHz is that between the elec¬ 
tron and the ^^N nuclear spin [b^ [t^- The magnetic 
fields Bx{t) and By(t) are time-dependent external con¬ 
trol fields. The random field B{t) represents the effect of 
the bath spins with correlation function [5l|, [7§-[8lj 

C{tf) = {B{t)B(t’)) 

= ( 6 ) 

where b that can be regarded as the field inhomogeneity 
in the effective semi-classical random noise model is as¬ 
sociated with the average coupling strength between the 
electron spin and the spin bath, and Tc is the correlation 
time of the spin bath. 

We will consider first the case where the quantum gate 
operations of the system qubits are influenced by the in¬ 
teractions with nearby noise qubits with a low number 
of degrees of freedom [bO, El, HI]- In this case, the most 
general approach is to describe the dynamics of both the 
system and noise qubits and their interactions in Hamil¬ 
tonian unitary approach and then perform the QOG cal¬ 
culations for the quantum gate operations. As the num¬ 
ber of degrees of freedom of the noise qubits increases, 
the computation of this approach becomes expensive and 
challenging. Our strategy is that we treat the spins which 
are close to and have large hyperfine interactions with 


the system qubit (s) as the noise qubits and include them 
and their interactions into the unitary Hamiltonian, and 
then take the many randomly distributed distant spins 
as a spin bath (environment) whose average effect on the 
dynamics of the system qubit (s) is obtained by tracing 
out the environment degrees of freedom (or more pre¬ 
cisely by performing an ensemble average over the clas¬ 
sical random noise) using the master equation approach, 
of the reduced density matrix. So our QOG treatment 
can simultaneously deal with the effects of leakage states, 
a few nearby noise qubits and a (spin) bath. 

III. QUANTUM MASTER EQUATION 

Since the ensemble average effect of the semi-classical 
random noise treatment of the spin bath in the absence 
of the external control fields can be modeled by a pure 
dephasing open system model [sJl, we employ the per¬ 
turbative time-local non-Markovian master equation to 
describe the coherent control and decoherence dynamics 
of the NV center electron qubit in our QOG study. Thus, 
following the the standard perturbation theory with the 
Born approximation, we can write the time-local non- 
Markovian master equation for the reduced system den¬ 
sity matrix as [SSl - IS^ 

^ = CsP{t) + [CzD{t)p{t) + h.c.] , (7) 

— t —t 

where Cg = — [Hs{t),»] and Cz = ^ [<S'z,*], and the 
dissipator D can be written as [s^ E3| 

D(t) = t')Ug{t, t')Szdt', (8) 

where the propagator superoperator Us(t,t') = 
T+e-^t' denoting the time-ordering 

operator. The symbol h.c. in Eq. © denotes the 
hermitian conjugate of its previous term. To solve the 
master equation 0 directly, one would need to evaluate 
the term lAs{t,t')Sz and then perform the integration 
for the dissipator D[t) of Eq. ([5]). This procedure is 
often numerically inefficient. Instead, since the bath 
correlation function given in Eq. (|6]) is in an exponential 
form, one can take the time derivative on Eq. ® 
and obtain straightly the differential equation for the 
dissipator D{t) as 

^^D{t) = -ih^Sz+{Cg{t)-^^D{t). ( 9 ) 

Equations o and 0 form a coupled set of inhomoge¬ 
neous differential equations, and one can use the Runge- 
Kutta method to solve these equations numerically. 

For the convenience of numerical computation, we 
transform the density matrix p into a column vector 
p and in this case Eq. © becomes p{t) = A{t)p{t), 
where A{t) is the corresponding operator associated with 



4 


Eq. 0 in column vector representation. It can be shown 
that the effective propagator U(t) defined by the relation 
p{t) = U{t)p(b) satisfies: 

U{t) = A{t)U{t) with C7(0) = If/, (10) 

where J\f is the dimension of U and Ij/ denotes the JVxAf 
identity matrix in the operator dimension of U. 


IV. GATE ERROR AND OPTIMAL CONTROL 
ALGORITHM 


where PU^ij) are ub xub matrix partitions of PU in the 
computational state basis of the system qubits, and Gij 
are the scalar matrix elements of the target operation G 
[S^ . Note that because the projected propagator PU is 
no longer unitary when the effect of the leakage states 
or/and the open system environment is considered, the 
second term in the error function m takes the form of 
Tr[{PU)'^PU]/{2N). This term approaches (1/2) for an 
ideal, closed composite system when PU —> C/ is unitary, 
and in this case the error function K (TT^ reduces to that 
of Ref. m. 


We describe below briefly the error function of the 
gates and the optimal control algorithm based on the 
Krotov method that we adopt for our calculations. The 
noise qubits are regarded as an effective small environ¬ 
ment interacting with the system qubits that serve as a 
register for quantum information processing. The imple¬ 
mentation of quantum gates in the presence of only a few 
noise qubits using the Hamiltonian unitar y a pproach for 
QOC calculations has been investigated [8^ |8^. Here 
besides a few noise qubits, the leakage state and the spin 
bath are also considered. We thus define the error func¬ 
tion K as the distance measure between the associated 
propagator PU{T) at the final time T of the compos¬ 
ite system and the target gate G in the column vector 
representation as follows [H, : 

K = Xn nun {II PU - G ® 4) |||| = 1^^} . (11) 

Here P denotes a projection operator to project the prop¬ 
agator U onto the composite subspace spanned by the 
tensor product of the system qubit computational basis 
states and the noise qubit basis states n_B is the 

dimension of the Hilbert space for the noise qubit sub¬ 
system, and $ is an arbitrary unitary acting only on the 
noise qubit Hilbert space. The symbol || • ||_f stands for 
the Frobenius matrix norm: || A ||f= (TrA^A)^/^, and 
Aat = ^ is a normalization factor to keep the value of 
K in the range [0,1] with N the dimension of PU. The 
squared norm of || PU — G 0 $ Hi. is minimized over 
the set of all possible unitary $ because we do not care 
what the evolution of the noise qubits is as long as the 
target gate G to be implemented can be achieved. The 
fidelity defined a.s F = 1 — K with K given in Eq. (HU is 
introduced to measure how well PU{T) approaches the 
target gate G at the final gate time T. The error function 
similar to K of Eq. (HU for a closed composite system 
consisting of the system qubits and a few noise qubits was 
defined and simplified to a computable form in Ref. |H. 
Here, we generalize the expression of the error function 
K for an open system with leakage states and/or a spin 
bath in a computable form as (cf. [ssjl l 

K=l + ^Tr[(PG)tpt/] - iReTrv^. (12) 
Here in Eq. HU 

ns 

Q = Y. G*^iPU\,,p (13) 


In realistic control problems, it is desirable that the 
optimal control sequence can provide the highest quality 
(fidelity) with minimum energy consumption. Therefore, 
we define the objective function to be maximized for our 
optimal control problem as 

j = F- £ (14) 

with a weighting function X{t) > 0 adjusted and chosen 
empirically. Here the reference field £{t) is chosen to be 
the control field of the previous iteration [H, . In this 

case, when the iterative procedure approaches the opti¬ 
mal solution, the change in the control field is minimal 
or vanishing. Therefore, this choice of the reference field 
£{t) ensures that the iterative method is found to increase 
the total objective J of Eq. HU by increasing the gate 
fidelity F = 1 — K rather than reducing the total control 
pulse energy. 

The iterative algorithm of the Krotov method used 
in our optimal control study of implementing quantum 
gates is described briefly as follows [sUIS^ . 1^ . 1^ . 1^ . 
pTfil . [ 7 ^ . (1) An initial guess for the values of the con¬ 
trol parameters e®(t) is randomly chosen [here the con¬ 
trol parameters £i{t) can be the externally applied ac 
magnetic fields Bi{t) with i = x,y components]. (2) 
The evolution propagator 17[e°(t)] is evolved forward in 
time until t = T using Eq. (HOI) . (3) An auxiliary func¬ 
tion S[e|(t)], / = 0 for the first iteration, is evolved 
backward in time until t = 0 using the equation of 
motion B{t) = —B{t)A{t) and the boundary condition 

dK dK 

B(T) = ~ ■ The explicit form of can be 

found in the Appendix (4) The updated propagator 
U[£i'^^{t)] is propagated again forward in time, while the 
control parameter £i{t) is updated iteratively with the 

rule£|+^(t) = e|(t)-k^Re|Tr[HJ (t)^|^C/J+i(t)]|. 

(5) Steps (3) and (4) are repeated until either the error 
is smaller than a preset value or the ratio of the im- 

proved error at the next iteration, -, is rather 

small. After a sufficient number of iterations, the algo¬ 
rithm converges and the fidelity F in the objective func¬ 
tion of Eq. HU reaches asymptotically a maximum value 
of -Fmax- 
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V. RESULTS AND DISCUSSION 
V.l. Optimal control for single-qubit gates 

The target quantum single-qubit gates considered in 
our investigation are ^-gate, A-gate, Hadamard gate {H 
gate), phase gate and tt/S gate on the NV center elec¬ 
tron spin. We restrict our investigation to have a control 
magnetic field only in the x-direction, i.e., Bx{t), for the 
implementation of these single-qubit gates. A static mag¬ 
netic field Bz = 500 G is applied to split the Ws = ±1 
states of an NV center electron spin; the rus = 0 state 
and iTis = — 1 state are chosen as the system qubit states, 
and the rus = 1 state is treated as a leakage state. We 
will investigate the case where there are two nearby noise 
qubits, the nuclear spin and nuclear spin, and 
the strength of the hyperfine interaction between the NV 
center electron spin and the nearest-neighbor nuclear 
spin is comparable to the Zeeman splitting of system 
qubit states. 

A quantum Z-gate in the absence of the noise qubits 
can be realized by free evolution of the system qubit with 
high fidelity (or error smaller than 10“® ). However, the 
Z-gate errors K obtained by free evolution taking the en¬ 
ergy shift by the and ^®N noise qubites into account 
at operation times of 0.34 ns, 7.45 ns and 15.2 ns are all 
greater than 6.0 x 10“^. On the other hand, using op¬ 
timal control method with an extra control field Bx{t) 
gives, for example, a Z-gate error K for an operation 
time of 10 ns (or 0.01/is) to be 6.0 x 10“'^, at least one 
order of magnitude better than that by the free evolution. 

For the implementation of an X gate, traditionally a tt 
pulse with frequency w in resonance to the energy split¬ 
ting between the two computational states of the system 
qubit (in our case w = 1.343 • 27r GHz at Bz = 500 G) 
is used to induce the qubit transition. As fast quan¬ 
tum gates are favorable for the purpose of quantum in¬ 
formation processing, we set the operation times of the 
single-qubit gates to be in the order of O.Ol^s. To be 
operated in such a short time, an X gate implemented 
by a TT pulse would require a certain pulse strength B^q 
of Bx{t) = Bxo cosujt. This in turn limits the fidelity of 
the gate as a relatively strong pulse strength B^o would 
cause possible transitions to the TOs = -1-1 leakage state. 
One can see this from Fig. Ha) that the gate errors of 
X gates implemented by tt pulses for different operation 
times in the absence of any noise qubit and decoherence 
represented by the green circles (guided by the green dot¬ 
ted line) are all larger than 4 x 10“^. In the presence of 
the nearby noise qubits, nonsecular parts of the hyper¬ 
fine couplings to the system qubit will cause an addi¬ 
tional error to the X gate implemented by a tt pulse. The 
purple, blue and red circles (dashed lines) in Fig. Ha) 
represent the errors of the X gates implemented by tt 
pulses at different operation times for the cases where 
the system qubit interacts with only the ^®N, only the 
^^G, and both the ^®N and the ^^G noise qubits, respec¬ 
tively. One can see that the gate errors of the blue and 
red dashed lines overlapping with each other are consid¬ 


erably larger than those of the purple and green dashed 
lines also overlapping with each other in the short gat¬ 
ing time regime but deviating a little bit at large gating 
times. This is because the hyperfine interaction of the 
^^N nuclear spin with the system qubit is about one to 
two orders of magnitude smaller than the control field 
strength or interaction for the X gates with operation 
times shown in Fig. Ha), and thus the presence of the 
^®N nuclear spin does not introduce substantial error in 
the X gate. On the other hand, the hyperfine interaction 
of the ^^G nuclear spin with the system qubit is compara¬ 
ble to the control field interaction and thus considerably 
larger gate error is introduced by ^^G nuclear spin than 
that by the leakage state in the ideal case or than that 
by the ^^N nuclear spin. 

In contrast, the QOG based on the Krotov optimiza¬ 
tion method enables us to achieve a high-fidelity X gate 
with significant improvement in gate error. As shown in 
Fig-Ha), the X gate errors obtained by the optimal con¬ 
trol method in the presence of both the ^^N and the ^^G 
noise qubits represented by the red squares (solid lines) 
are two orders of magnitude smaller than the gate error 
obtained by directly applying tt pulses even in the ideal 
case (green dotted line). Figure Hb) shows a typical opti¬ 
mal control field sequence for a X gate with an operation 
time of 0.01 /iS. We have also performed optimal control 
calculations for the single qubit gates of the Hadamard 
gate {H gate), the phase gate and the tt/S gate in the 
universal quantum gate set for FTQG. The errors of the 
single-qubit gates with operation times all set to 0.01 fis 
are summarized in Table [T] Because the gate times of the 
single-qubit gates are all set to a relatively short time of 
0.01 IIS, we find that the calculated optimal control pulse 
sequences are robust (i.e., the gate errors do not increase 
appreciably) in the presence of a spin bath with a wide 
range of realistic experimental parameters for the bath 
correlation function defined in ([B]) [13, HH, Ull ■ The effect 
of a spin bath will be explicitly discussed for two-qubit 
CNOT gates that have longer gate times. 


V.2. Optimal control for CNOT gates 

We next describe the implementation of the two-qubit 
GNOT gate in the discrete set of universal gates. We 
choose the electron spin of the NV center as the con¬ 
trol qubit and the ^^G nuclear spin as the target qubit. 
The ^^N nuclear spin associated with the NV center acts 
as a noise qubit influencing the electron qubit. Besides, 
the bath of other distant spin impurities causing the de¬ 
coherence to the electron qubit is taken into account. 
A schematic illustration of the whole system we con¬ 
sider is shown in FigHb). The experiments carried out 
on a single electron and a single ^^G nuclear spin of a 
NV center for the implementation of a two-qubit condi¬ 
tional rotation (GROT) gate realized by a conditional 
radio-frequency tt pulse have been reported [1, . This 

GROT gate combined with single-qubit 2 rotations can 
perform a CNOT gate up to a global phase factor. De- 
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(a) 



(b) 


FIG. 2. (Color online) (a) X gate errors versus operation 
times at a static magnetic field Bz = 500 G for different 
noise-qubit scenarios and operation schemes. The colors of 
purple, blue and red represent the cases where the system 
qubit (electron spin) interacts with only the ^®N, only the 
and both the and the noise qubits, respectively. 
The open circles with dashed lines represent X gates imple¬ 
mented by conventional tt pulse, while the solid squares with 
solid lines denote implementation by optimal control pulse se¬ 
quence. The open green circles with a dashed line stand for 
the X gate error of an ideal NV electron spin without any 
noise qubit nearby, (b) A typical Optimal control pulse se¬ 
quence of Bx field for a gate operation time of 0.01 /rs at a 
static magnetic field of Bz = 500 G. 


composing a CNOT gate or a general gate operation into 
several single-qubit and entangled two-qubit operations 
in series makes its operation time normally longer and 
its overall gate error larger as the gate errors of the de¬ 
composed gates will accumulate. In contrast, the optimal 
control method has the great advantage of enabling the 
implementation of a CNOT gate or other general quan¬ 
tum gates in a single run of pulse or in a single pulse 
sequence by simply setting the target operation to be the 
CNOT gate or the general quantum gate one wishes to 
implement [fill [zij . Simulating this 7r-pulse approach of 


TABLE I. Summary of the QOC gate errors. The bath corre¬ 
lation time for the CNOT gates is Tc = 25 /rs, and the aver¬ 
age system-bath coupling strengths (or field inhomogeneities) 
are b = 38.46 /rs”^ for the gate time T = 0.125 /rs and 
6 = 80 /is“^ for the gate time T = 0.05 /rs, respectively. 


gate type 

noise qubit (s) gate time T gate error K 

Z gate 

«N 

15n,13c 

0.01 fis 
0.01 /is 

8.0 X 10”'’ 

6.0 X lO"** 

X gate 

15n,13c 

0.01 /is 
0.01 /is 

5.5 X 10”® 

3.9 X 10"® 

H gate 

15n,13c 

0.01 /is 
0.01 /is 

2.4 X 10”® 

1.0 X lO"** 

Phase gate 

15n,13c 

0.01 /is 
0.01 /is 

2.5 X 10"® 

4.3 X 10”® 

n/8 gate 

15n,13c 

0.01 /is 
0.01 /is 

2.3 X 10"® 

1.4 X 10”'^ 

CNOT gate 

■tsn 

i®N, bath 

0.125 /rs 
0.125 fjs 

1.7 X 10"'‘ 

6.5 X 10"^ 

CNOT gate 

"f®N 

i®N, bath 

0.05 /is 
0.05 /is 

4.3 X 10-'* 

4.3 X 10"'^ 


implementing the CROT gate using Bx{t) = Bxocosojt 
and By (t) = Byo sin ujt with the field strengths Bxq and 
RyO) we evaluate the CROT gate error K and compare 
it with our result of CNOT gate by the optimal control 
method. 

We discuss the implementation of a CNOT or CROT 
gate for three different cases. The first case is the ideal 
case without any noise, the second case includes the ef¬ 
fect of the ^^N noise qubit, and the third case considers 
both the ^®N noise qubit and a spin bath. Here applying 
a static magnetic field of Bz = 1000 G, results in a rel¬ 
atively large energy splitting between iris = 1 (leakage) 
state and m = 0, —1 (computational basis) states as com¬ 
pared to the energy splitting between the computational 
basis states of the TOs = 0 and TOs = — 1 states for the NV 
electron spin. At first, the NV electron spin is regarded 
as an ideal two-level qubit system for the gate implemen¬ 
tation, and then we will treat the NV electron spin as a 
spin-1 three-level system. This enables us to exam how 
well the optimal control field sequence obtained for the 
two-level case performs in the more realistic three-level 
case (i.e., to see the leakage effect). We then employ the 
QOC theory to find new control sequence for the three- 
level NV center electron spin taking the leakage state of 
TOs = 1 into account in order to reduce the gate error or 
infidelity. 

Ideal case: Next we show that the optimal control 
theory can achieve a CNOT gate with a better fidelity 
than simply applying a conditional tt pulse (CROT). The 
CNOT or CROT gate errors as a function of the gate op¬ 
eration time for different cases are illustrated in Fig. |3l 
The solid lines and dot-dashed lines represent the gate 
errors obtained for the ideal case and the case includ¬ 
ing noise qubit ^®N, respectively. The red lines represent 
the CROT gate errors obtained by applying a condi- 
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FIG. 3. (Color online) CNOT or GROT gate errors versus 
operation times at a static magnetic field = 1000 G for 
different operation schemes: Conditional tt pulse scheme (in 
red) and optimal control scheme taking the NV center electron 
spin as a two-level system (in blue) and three-level system (in 
purple). The solid lines and dot-dashed lines represent the 
gate errors obtained for the ideal case and the case including 
noise qubit respectively. 


tional TT pulse taking iris = 0,-1 states of the NV center 
as the system qubit states. The blue lines and purple 
lines represent the CNOT gate errors obtained by using 
the optimal control theory taking the NV center elec¬ 
tron spin as a two-level system and three-level system, 
respectively. Note that different from the gate error of 
the CROT implemented by a conditional tt pulse, the 
QOC gate errors are calculated by substituting the pulse 
sequence obtained in respective methods into the Hamil¬ 
tonian of the three-level NV center system including the 
^®N noise qubit (but not including the bath spins) for the 
propagator U in Eq. For the idea case (solid lines), 
when the gate operation time becomes shorter, the gate 
error using a tt pulse becomes larger. The reason is that 
the requirement for a shorter duration of a tt pulse means 
a stronger field strength, which in turn induces a larger 
transition probability to the leakage state. In contrast, 
the QOC theory can do a much better job when the op¬ 
eration time is short as illustrated in Fig 121 However, as 
the operation time increases, the error for the optimal 
control case of treating the electron spin as an effective 
two-level system becomes larger, even larger than that of 
using a tt pulse (see Figl21). Therefore, employing QOC 
for the three-level electron spin system is necessary. It is 
quite obvious from Fig l3] that employing the QOC theory 
by treating the electron spin as a three-level system (pur¬ 
ple solid lines) gives gate errors one order to two orders 
of magnitude lower than those of the two-level case (blue 
solid line). 

Effect of a noise qubit: In the case of including the 
^®N noise qubit, driving the system by a tt pulse does not 
work well anymore (see the red dot-dashed line). How¬ 
ever, taking the noise qubit into the QOC optimization 



(G) 

(a) 




(b) 


FIG. 4. (Golor online) (a) Optimal control CNOT gate errors 
for different values of the static magnetic field Bz taking the 
NV center electron spin as a two-level system (blue solid line) 
and three-level system (red dot-dashed line). The operation 
time is 0.05 fis. (b)Optimal control pulse sequences of B^it) 
and By{t) for Bz = 500 G (blue solid lines) and Bz = 1000 G 
(red dot-dashed lines) 


consideration lowers the error or infidelity a little bit as 
compared to the case without doing so. One can see 
this from Figl3]that the blue and purple dot-dashed lines 
are slightly lower than their corresponding blue and pur¬ 
ple solid lines, respectively. The reason is that the ^^N 
noise qubit here also serves as an ancilla qubit that allows 
some probability to get out of the computational space 
temporarily but return to the computational space with 
higher fidelity at the end of the operation for the optimal 
control CNOT gate. Again, treating the electron spin 
as a three-level system is a more effective strategy to per¬ 
form CNOT gate and it gives an error K Ri 1.7 x 10“^ 
in the presence of the ^^N noise qubit for a gate time 
of 0.125 ^s. Figure SKa) investigates the optimal con¬ 
trol CNOT gate error as a function of the strength of 
the external static magnetic field for a gate opera¬ 
tion time of 0.05 fis. The blue solid line and red dot- 
dashed line represent the cases of treating the electron 
spin as a two-level system and a three-level system, re¬ 
spectively. When the static magnetic field decreases. 























FIG. 5. (Color online) CNOT gate error verse the average 
system-bath coupling strength (field inhomogeneity) h for gat¬ 
ing time of = 0.05 ^s (blue solid line) and 0.125 ^s (red solid 
line). The corresponding dashed lines represent the gate er¬ 
rors calculated by taking the OCT pulse sequences obtained 
for the case without the spin bath (i.e., setting 6 = 0) as our 
initial guess for optimal control in the open system. 


the gate error represented by the blue solid line becomes 
large. This is because at a lower magnetic field Bz the 
energy separation between the qubit state mg = — 1 and 
the leakage state Wg = 1 become smaller, and thus the 
approximation of treating the NV electron spin as a two- 
level system of Wg = 0, —1 is not very good. However, 
the optimal control of the three-level case still works with 
error K k, Qx 10“"^ even when the magnetic field is as low 
as Bz = 500 G. The optimal control field sequences by 
treating the NV electron spin as a three-level system for 
Bz = lOOOG (red dot-dashed lines) and Bz = 500G (blue 
solid lines) with the operation time of 0.05 p,s are shown 
in Fig. SKb). At Bz = 500 G {Bz = 1000 G), the energy 
separation between the electron qubit states of rUg = — 1 
and rrig = 0 is much larger than (comparable to) the 
energy splitting between the states | 1 ,),) and |l,t) (due 
to the hyperfine interaction), where |),) and |t) represent 
the nuclear spin states, and |1) and |0) represent 
the electron qubit states. As a result, there are appar¬ 
ently two different oscillating components for the case of 
Bz = 500G (blue solid lines) as compared to the case of 
Bz — 1000 G (red dot-dashed lines) in Fig. SKb). The 
major frequency of the fast oscillating component comes 
from the high frequency (energy separation) of the elec¬ 
tron qubit, while that of the slow one matches the energy 
splitting between states |1,),) and |l,t)- 

Effect of a spin bath: Next, we take into account not 
only the noise qubit but also the bath of the distant spin 
impurities, which is modeled as an effective classical ran¬ 
dom field acting on the NV electron spin with correlation 
function given by Eq. The solid lines in Fig. [S] show 
the dependence of the CNOT gate error on the bath pa¬ 
rameter b for the correlation time Tc = 25 /iS. The param¬ 
eter b defined in the bath correlation function, Eq. ®, 
is related to the average system-bath coupling strength 



FIG. 6. (Color online) CNOT gate error verse the correlation 
time Tc for the case of 6 = 38.46 /is“^ and the gating time 
T — 0.125 fis and the case of 6 = 80 fis~^ and the gating 
time T = 0.125 fis) The dashed line represents gate error 
calculated by taking the OCT pulse sequences obtained for 
the case without the spin bath (i.e., setting 6 = 0) as our 
initial guess for optimal control in the open system. 


or field inhomogeneity in the random field model. The 
dashed lines in Fig.[S]represent the gate error obtained by 
applying the optimal control sequences obtained for the 
case without the spin bath to the master equation or the 
effective evolution equation ([T0|) in the presence of the 
bath. This allows us to study how the gate errors due to 
the spin bath are improved in the open system by QOC. 
One can see from Fig. [S] that for Tc = 25 fiS unless the 
parameter b is more than one order of magnitude larger 
than the typical value of 6 = 3.846 /is“^ extracted from 
the experiments [SlEllIl , the influence of the bath 
on the gate with a short operation time of 0.05 ps (blue 
dashed line) is not significant. For the longer gate oper¬ 
ation time of 0.125 ps red dashed line), the influence of 
a bath becomes appreciable when 6 > 5 ps~^, and as the 
coupling strength b become stronger, the error increases. 
However, the optimal control CNOT gate with operation 
time 0.05 ps (blue solid line) can sustain much larger val¬ 
ues of the field inhomogeneity up to 6 = 80 ps~^ and the 
gate error can be maintained almost the same as if the 
spin bath were not present. 

The relation between the gate error and the bath 
correlation time Tc for the average coupling strengths 
of 6 = 38.46 ps~^ and the operation time of 0.125 ^s 
(6 = 80 ps~^ and operation time 0.125 ps) is shown in 
red solid line (blue solid line) in Fig. [0] To investigate 
how much QOC theory improves the gate fidelity in the 
NV center system with a spin bath, we take the opti¬ 
mal pulse sequences in the absent of the spin bath as 
our initial guess for the control fields, and the gate errors 
before the optimal control iterations are shown in dot- 
dashed lines in FiglHl The gate errors (solid lines) after 
the iterations are terminated are larger for small values 
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of Tc, and decrease as Tc increases. The gate errors (solid 
lines) stay nearly constant as Tc > 5 fj,s. When the bath 
correlation time Tc is small, the memory effect of the bath 
is weak and the bath is close to a Markovian bath. In this 
case, it is hard to revise the bath contribution since the 
decay rate approaches to a steady value in a very short 
time. In contrast, when Tq is large, due to the relatively 
long memory effect of the bath, the optimal control fields 
are able to counteract the influence from the bath. 


VI. CONCLUSION 


We have found the control sequences of fast and high- 
fidelity single-qubit and two-qubit quantum gates for 
electron and nuclear spins of a NV center in diamond 
using the QOC theory. A CNOT gate or other gen¬ 
eral quantum gates operation can be implemented in a 
single run of pulse sequence using the optimal control ap¬ 
proach rather than being decomposed into some entan¬ 
gled two-qubit and several single-qubit operations in se¬ 
ries by composite pulse sequences. The (non-Markovian) 
external environment effects on system qubits are parti¬ 
tioned into two kinds, one from a few nearby noise qubits 
and the other from a bath of distant spins. Table U sum¬ 
marizes the gate time and gate errors calculated with 
realistic experimental parameters for the cases consider¬ 
ing the effect of the noise qubits, the leakage state and 
the effect of a spin bath. These gate errors are below 
the recent model of error threshold 10~^ 17^ (10“^ if 
surface code error correction is used |73l475l| l required 
for FTQC. One can estimate the logical error rate that 
our gate error K corresponds to for a given type of er¬ 
ror correction of FTQC. Suppose that one implements 
a logical qubit of the surface code error correction on a 
two-dimensional array of physical qubits with array size 
(distance) d [corresponding to the number of physical 
qubits being Ug = {2d — 1)^] as in Ref. [t^. The num¬ 
ber of physical qubits or the array size d needed to de¬ 
fine a logical qubit to meet a required logical error rate 
is dependent strongly on the error rate in the physical 
qubits. By taking the error probability per step to be 
the worse gate error K of 6.0 x 10“"^ in Table U which is 
smaller than the error threshold rate of pth = 0.57% of 
the surface code, the estimated logical error rates using 
Fig. 4 or more precisely Eq. (11) of Ref. [t^ are about 
3.5 X 10“^ for array size d = 5 and about 3.7 x 10“® for 
array size d = 7. When the per-step error rate is smaller 
than the error threshold rate pth, the logical error rate 
falls exponentially with the array size d. Thus one ob¬ 
tains the logical error rate for array size d = 25 to be 
about 5.8 X 10“^^, sufficient to perform Shor’s algorithm 
for factoring a 2000-bit number into its primes with a 
reasonable chance of success [^ . 

There seems to be challenges for the experimental im¬ 
plementations and applications of the QOC theory, such 
as imprecise knowledge of the quantum system’s param¬ 
eters and how to generate the complex optimal con¬ 
trol pulses in reality. Fortunately, commercial devices 


for generating arbitrary wave forms or complex signals 
in a time scale of sub-nanoseconds to nanoseconds are 
available now and may solve the challenge of generating 
complex pulse sequences. Besides, a hybrid open-loop- 
closed-loop optimal control method called adaptation by 
hybrid optimal control (Ad-HOC) method [S^ designed 
to overcome not only the problem of inaccurate knowl¬ 
edge of the system parameters but also shortcomings of 
the assumed physical model and errors on the control 
fields has been recently proposed. The closed-loop pulse 
calibration of Ad-HOC, similar to adaptive model-free 
feedback control (also referred to as closed-loop labora¬ 
tory control or learning control) Bin, uses the physi¬ 
cal system itself as a feedback to calibrate control pulses 
and optimize their performance. Two similar closed-loop 
methods for optimizing quantum control in experimen¬ 
tal systems have also been put forward recently: the 
method of optimized randomized benchmarking for im¬ 
mediate tune-up (ORBIT) and the method of adap¬ 
tive control via randomized o ptim ization nearly yielding 
maximization (ACRONYM) [^. In principle, the sys¬ 
tem response and control pulses can be calibrated and 
improved using closed-loop optimization where measure¬ 
ment data are efficiently obtained with Nelder-Mead al¬ 
gorithm [^, or stochastic optimization algorithm 
[93 and subsequently fed back to the system op¬ 

timizer to improve the pulses without precise knowledge 
of the system. These developments make the QOC the¬ 
ory practical and useful to construct the initial pulse se¬ 
quences for experimentally closed-loop optimization [s^ . 
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Appendix A: Explicit form of dK/dU 


Employing the algorithm of the Krotov optimization 
method, we evolve backward in time the auxiliary func¬ 
tion B{t) with the boundary condition B{T) = 
Therefore, we need to find the explicit form of The 
procedure and calculation presented below follows those 
in Ref. [H 

In our case, the error function K is defined through 
Eqs. (ini) and (131). The matrix elements of can be 
written as 


f dK\ _ dK 

\dPU) dPUab ^ ^ 

where the matrix indices a and b range form 1 to V 
and PUab is the matrix element (a complex scalar vari¬ 
able) of PU. Let X — PUab , and Z{x) = 
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be a matrix function of the variable a;, and y{Z(x)) = 
Tr-yQlQ = TyZ{x) be a scalar function. Thus, the 
error function (HU in terms of the new variables is 
K = i + ^Tr[{PU)^PU] — -^y. Note that \/Q^Q is not 
an analytic function of PUab but can be expressed as an 
analytic function of PUab and PU*^^. Consequently, when 
one differentiates \/Q}Q with respect to PUab, PU*i^ and 
subsequently are treated as constants. Thus 


dK 

dx 


^PUL 

^PUL 

^PUL 

^PUL 

-^PPba 


N dx 

'E 


N 


k,k' 


-E 

N ^ 


k,k' 


-Tr 


N' 


dy dZk,k' 
dZk,k' dx 

dy dZl,k 
dZk,k' dx 

' dy dZ"’" \ 

^dZ dx J 

'dZ'^'S 

. dx J ' 


(A2) 


Here the second line in Eq. (IA2I) follows from the use of 
the chain rule, and the last line follows from the property 


of — = - -^ = I. Then from the definition of Z 

dZ dZ 
and x, we find that 


dZ _ d^jQ'^Q _ 1 

dx ~ dPUab ~ 2^^ ^ dPUab ■ 


(A3) 


Since Tr 
into Eq. 



■ = Tr—, inserting the trace of Eq. (IA3I) 
dx 

, we finally obtain the explicit form of 


\dPuJ^ 


2N 


PUL - Tr 




(A4) 


where 


dQ 

dPUab 


= |a mod ub) {b mod ns| 


(A5) 


can be simply obtained by the definition of Q in Eq. m- 
Here ia and jb denote the smallest integers greater than 
or equal to a/ns and b/uB, respectively, and Gi^j^ are 
the matrix elements of the target gate. In Eq. (lAhl) . 
the states \a mod ub) {b mod nsj are the elements of the 
orthonormal basis matrix of the noise qubit (s). 
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